Helical superconducting black holes 

Aristomenis Donos^ and Jerome P. Gauntlett^ 
'Theoretical Physics Group, Blackett Laboratory, Imperial College, London SW7 2AZ, U.K. 

We construct novel static, asymptotically AdSf, black hole solutions with Bianchi VIIo symmetry 
that are holographically dual to superconducting phases in four spacetime dimensions with a helical 
p-wave order. We calculate the precise temperature dependence of the pitch of the helical order. At 
zero temperature the black holes have vanishing entropy and approach domain wall solutions that 
reveal homogenous, non-isotropic dual ground states with emergent scaling symmetry. 
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INTRODUCTION 

The AdS/CFT correspondence is a powerful tool to 
analyse strongly coupled quantum field theories and there 
has been a surge of activity aimed at finding possible 
applications both to condensed matter systems and to 
QCD. One focus has been to holographically realise var- 
ious kinds of phases via the construction of fascinating 
new classes of black hole solutions, which are also of in- 
terest in their own right. 

An important development was the discovery of black 
brane solutions that are holographically dual to super- 
conducting phases, or more precisely, superfluid phases 
[1] . These black holes carry a halo of charged hair which 
spontaneously breaks a global abelian symmetry of the 
dual field theory. In the simplest examples, the charged 
hair is in the guise of a bulk scalar field, correspond- 
ing to a scalar order parameter in the dual field theory, 
and hence an s-wave superconducting phase, p-wave su- 
perconducting phases, in which the order parameter has 
angular momentum I = \, have also been realised [2ll3j. 

Spatially modulated phases are also widely seen in Na- 
ture. The order parameters for these phases are associ- 
ated with non-zero momentum and spontaneously break 
some or all of the translation invariance. Common ex- 
amples in condensed matter include spin density waves 
and charge density waves, while QCD at high baryonic 
density is anticiapted to be in a chiral-wave state [S]. 
Spatially modulated phases that are also superconduct- 
ing are possible [1] and such FFLO phases have been 
argued to be realised in several systems [6j . Of particu- 
lar interest here are p-wave superconducting phases with 
a helical order. In these phases the I = 1 order parameter 
points in a given direction in a plane which then rotates 
as one moves along the direction orthogonal to the plane. 
They have been discussed, for example, in the context of 
non-centrosymmetric heavy fermion compounds [7]. 

Holographic studies of spatially modulated phases were 
initiated in [5]. The purpose of this letter is to present 
the very first construction of fully back-reacted black hole 
solutions that are holographically dual to spatially mod- 
ulated phases which are, moreover, superconducting. We 
will consider a class of gravitational models in Z) = 5 that 
couple a metric with a gauge-field and a two-form poten- 
tial, which have been shown, using a linearised analysis. 



to admit black brane solutions which are dual to p-wave 
superfluid phases with a helical order in d = 4 [Sj . 

Our construction of the new black holes allows us to 
show that the helical p-wave superconducting phase is 
thermodynamically preferred and that the phase transi- 
tion is generically second order. The helical order is fixed 
by wave-number fc, or equivalently a pitch p — 2n/k. We 
calculate k{T) and find it monotonically decreases down 
to a finite value as T — >■ 0. We find that the solutions have 
vanishing entropy density as T — )• 0. In this limit they 
approach smooth domain wall solutions, which we also 
construct, that interpolate between AdS^ in the UV and 
a new spatially homogeneous but non-isotropic ground 
state in the IR, of a type recently been discussed in TIT . 

There are other contexts in which spatially modulated 
black holes should exist but, in general, the construction 
will require solving non-linear partial differential equa- 
tions. By contrast, a key point here is that the black holes 
for the helical p-wave superconductors can be obtained 
by solving ordinary differential equations since they are 
static and also have a Bianchi VIIq symmetry. 



THE D=5 MODEL 

As in W we consider a D — 5 model coupling a metric 
to a gauge field A and a complex two-form C with action 



5 = 



d^a 



-g[R+ 12 -\F^,F^^-~\C,,,C^''' 

+ 'ibri'^^''^'^ CtivHp„s\ , (1) 



where a bar denotes complex conjugation and the field 
strengths, using a form notation, are given by 



F = dA, 



H = dC + ieAAC. 



(2) 



This simple class of models, specified by the parame- 
ters m and e, is rather natural. The equations of motion 
admit a unit radius AdS^ solution with A = C = 0, 
which is dual to some putative conformal field theory 
(CFT). The massless gauge-field A is dual to a current 
in the CFT, corresponding to a global abelian symmetry, 
with scaling dimension A = 3. The two-form C sat- 
isfies a first-order equation of motion and is dual to a 
self-dual rank two tensor operator with A = 2 -f- |?7i|. In 



particular, this charged operator has I = 1 and thus pro- 
vides an order parameter for p-wave superconductivity. 
Such two-forms are common in Kaluza-Klein reductions 
from D=10 or D=ll supergravity. For example, when 
e — 1/VS, m = 1 precisely this model can be obtained as 
a consistent Kaluza-Klein truncation of type IIB super- 
gravity on S^ and moreover the operator in iV = 4 SYM 
dual to C is known j3l. 

We will study the CFT dual to the AdS^ vacuum at fi- 
nite temperature T and chemical potential fi with respect 
to the global abelian symmetry by constructing electri- 
cally charged asymptotically AdS^ black branes. 



BLACK HOLE SOLUTIONS 

The ansatz for all the black hole solutions that we con- 
sider is given by 



ds^ = -g f dt^ + g^^dr^ + h^Lol + r^ (e^" w^ -f e 
C = {icidt + C2dr) A 0^2 + C3 oji A wa , 
A — adt , 

where the one-forms lo^ are given by 

wi = dxi , 

W2 = cos {kxi) dx2 — sin (fcxi) dxy, , 

W3 = sin {kxi) dx2 + cos (kxi) dxs , 



'^'^uD 



(3) 



(4) 



and /, g, h, a, Ci and a are all functions of the radial 
coordinate, r, only and A; is a constant. Observe that k 
can be scaled out of the ansatz by scaling h but has been 
included for later convenience. The unit radius AdS^ 



vacuum solution can be obtained by setting g 



I 



1, h = r, a = a = Ci = 0. Notice that the ansatz (Is]) 
is static and in addition the constant t and r slices are 
spatially homogenous of, generically, Bianchi type VIIq. 
Substituting this ansatz into the D = 5 equations of 
motion we find that we can solve for ci and C2: 



„2q 



Cl 



C2 = 



1 






e "keacs + mhfgc^ 



fge^°'B + m?K^ 



ihc-i 



'kfgc'. 



(5) 



The remaining equations can be obtained from a one- 
dimensional action obtained by substituting the ansatz 
into the action associated with Q. 

The equations of motion admit the electrically charged 
AdS-RN black brane solution with a = Ci = 0, / = 1, 
h = r and 



M2 
3 



M 



(6) 



This solution approaches the unit radius AdS^ solution 
as r — >■ 00. The event horizon is located at r = r+ and 



the temperature is T = {Qr\ — /x^)/67rr_|_. This solution 
is dual to the CFT at finite chemical potential ^ and 
high temperatures. Clearly this phase is spatially homo- 
geneous and isotropic. It was shown in [S] that below a 
critical temperature, depending on the parameters tti, e, 
and the scale set by /i, this black hole is unstable to the 
formation of black holes that are dual to p-wave super- 
conductors with helical order. 

Let us first discuss the boundary conditions to be im- 
posed for the new black brane solutions. Regularity at 
the horizon demands that g{rj^) = a{rj^) = 0. We then 
find the solution at the horizon is specified by six param- 
eters: r+, /(r+), ft.(r+), a(r+), a'{r+) and C3(r+). As 
r — >■ cx) we approach AdS^ with asymptotic expansion 

5 = r2 (1 - Mr-^ + •••), / = /o (1 - Chr-^ + ■■■), 



h = r [1 + ciiT +•••), a = Car + ■ 
a = fo {^J. + qr^^ H ) , C3 = c„ 



"l™l 



(7) 



which is specified by eight parameters 
M, fo,Ch,Ca,fJ,,q,Cv and k. A number of comments 
are in order. Firstly, the fact that h ^ r implies that the 
wave-number k can no longer be scaled away. Secondly, 
the fall-off of C3 is chosen so that the charged operator 
dual to the two-form C has no deformation but can 
acquire, spontaneously ^ an expectation value propor- 
tional to c„ and spatially modulated in the xi direction 
with period 2-K/k. The holographic interpretation of 
the other UV parameters will be given below. Thirdly, 
there are two scaling symmetries of the differential 
equations which allow us to set ^ = fo ~ 1, and we will 
do so later (it is helpful to have them to discuss the 
thermodynamics). Finally, we have four second order 
differential equations for h, a, a, C3 and two first order 
equations for g, f and hence a solution is specified by 
ten integrations constants. On the other hand we have 
fourteen parameters in the boundary conditions minus 
two for the scaling symmetries. We thus expect a two 
parameter family of black hole solutions that can be 
specified by temperature T and wave number k. 



Action and Thermodynamics 

We analytically continue by setting t = —it and defin- 
ing / = —iS. The total action, including relevant 
boundary-terms, is given by 



^Tot 



I 



drd^x. 



-2K + & + ...]. 



(8) 



with goo = linir 



.9 



1/2 



fhr^ and K is the trace of the 



extrinsic curvature of the boundary at r — ?► 00. In ([8|) 
the ellipses refer to terms that will not contribute to the 
class of solutions that we are considering. The period 
of our Euclidean time is taken to be At and the tem- 
perature is then given by T = (/qAt)"^. Regularity at 



45'(47r)- 



We next define 



r = tj^ implies that T — -2^. 
the thermodynamic potential W — T[lTot\os = WV0I3, 
where [lTot\os is the on-shell action. Following the cal- 
culation in |llj we obtain the two equivalent expressions 

w = -M = e + 2/ig - Ts , (9) 

where we defined the entropy density s — Anr'^h-f-, s ~ 
3M+8ch and have set /o = 1- By calculating the on-shell 
variation of Ixot, including variations of /o, we deduce 
that w = 'w{T, /i) and the first law 



5w 



-s5T + 2q5^ . 



(10) 



Using ^ we also have 5e — T6s — 2fiSq. The identifica- 
tion of e with energy density is confirmed by computing 
the boundary stress-energy tensor [T^], again with /g = 1, 



3Af + 8ch , T, 



M + 8ch 



J-X2X2 ^^'^ 

Tx3Xa = M 



SCa COS {2kxi) 
8ca COS {2kxi) . 
8 Cq sin {2kxi) . 



(11) 



Observe that when Ca ^ the pressures in the X2, X3 
plane are spatially modulated as one moves along the Xi 
direction as one expects for helical order. Defining the 
average hydrostatic pressure, p, as minus the average of 
the trace of the spatial components we get p — M+8ch/^- 
The thermodynamically preferred black hole solutions 
that we construct will have Ca ^ and c^ — 0. Using ^ 
we conclude that this class satisfy the thermodynamic 
relation e + p = Ts — 2^q. 



Helical superconducting black holes 

The AdS-RN black brane solution ^ is unstable when 
e^ > m^/2 [9 . We will now consider the specific model 
with m — 1.7, e = 1.88 (for reasons we explain below) 
and set /x = /o = 1. For high temperatures we only find 
the AdS-RN black hole solution. The first new black hole 
solution appears at Tc ~ 0.0265 and for k — kc ~ 0.550. 
Holding this value of k fixed we construct these black 
hole solutions, numerically, all the way down to very low 
temperatures. Below Tc, as expected from the linearised 
analysis of [S] , there is a continuum of black hole solutions 
that appear with different values of k. By again holding 
k fixed we can construct each of these black holes too, 
down to low-temperatures. 

In fig. 1 we have summarised this new two-parameter 
family of solutions as well as displaying their free-energy 
w. All of these solutions have smaller free energy than the 
AdS-RN black brane solutions at the same temperature. 
At a given temperature T < T^. there is a one parameter 
family of black hole solutions specified by k and the one 
with the smallest free energy is depicted by a point on 




FIG. 1: The two-parameter family of helical 

superconducting black holes and their free-energy w. 

The red line denotes the thermodynamically preferred 

locus. The blue line is the free-energy of some domain 

wall solutions. The black lines are lines of constant T. 
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FIG. 2: Plots of c^ and wave-number k, which together 

fix the helical superconducting order, versus T for the 

thermodynamically preferred black hole solutions on the 

red line in fig. 1. The blue dots depict the quantities for 

the domain wall solutions. Note the scaled axes. 



the red line in fig. 1. Thus the one-parameter family of 
solutions specified by the red line characterises the ther- 
modynamically preferred solutions. Notice that we have 
a second order phase transition at T = Tc, k = kc and 
that as the temperature is lowered, the system smoothly 
moves between black hole solutions with different values 
of fc, all the way down to very low temperature where 
fc = fco ~ 0.256. In particular, the T = ground state 
remains spatially modulated. 

Interestingly, while the general two-parameter family 
of solutions have Ch 7^ 0, the solutions on the red line have 
(up to numerical accuracy) Ch = 0. In fig. 2 we have plot- 
ted, for the red line of solutions, the behaviour of Cy and 
wave-number fc, which, together, characterise the helical 
superconducting order, versus T. Near Tc we find the 
mean field behaviour Cy w 1.7 x lO^T^'^ (1 - T/Tcf^. 



New ground states at T = 

We are particularly interested in the T = limit of the 
thermodynamically preferred black hole solutions (the 
red line in fig. 1), which has fc = fco ~ 0.256. It is help- 
ful to first consider the T = limit of the whole class 
of black holes for general values of k. We find that they 
all approach a smooth domain wall solution that inter- 



polates between AdS^ in the UV and a new fixed point 
in the IR with a scaUng symmetry of a type that is very 
similar to those of [10] . Indeed we checked this exphcitly 
for the range 0.253 < k < 0.75 (going to smaller values 
of k becomes increasingly difficult numerically). 

To obtain this new fixed point solution, in ^ we put 

C3 = kcor, 



a = aor^ , 



h — kho , a — ao ■ 



(12) 



where L, ho^ao,ao, cq and z are all constant. By scaling 
t and xi we can set fo — k — 1. Notice that this ansatz 
corresponds to a solution invariant under the anisotropic 
scaling r — > X^^r, t — > A^i, 2:2,3 — > Aa;2,3 and xi — > xi. 
After substituting into the equations of motion we obtain 
a system of algebraic equations which can be solved. For 
the specific case of ttt- = 1.7, e = 1.88 we find 



z w 1.65..., 
ao« -0.380. 



L « 0.995 . . . 
ao« 0.265. 



hQ « 0.993 . . . , 
Co « 3.69... (13) 



We next construct domain wall solutions that interpo- 
late between this fixed point in the IR and AdS^ in the 
UV|15]. In the UV we continue to demand the expan- 
sion given by (It]). To obtain the IR expansion, we first 
consider perturbations about ( 12 ) of the form 



g = r^ (L + Xw^r') , / = /o r^-^ (l + Xw2r') , 

h = k [ho + Xw^r ) , a — uq + Xw4r , 

a = fodar'- (l + Xw^r^) , C3 ^ kc^r (l + Xw^r^) . (14) 

After expanding the equations of motion at first order in 
A we obtain a homogeneous linear system of equations 
E ■ w = where E is a 6 x 6 matrix that depends on 5. 
Demanding non-trivial solutions for w we determine the 
values of 5 by solving the polynomial equation |E| = 0. 
For the special case m = 1.7, e = 1.88 the modes with 
non-negative real parts have Sq = 0, (5i ~ 0.394, 82 ~ 
0.826, (53 « 0.847 and 5^ « 2.289 [H]. The IR expansion 
is then specified by a constant for each of these modes. 
Notice that the mode with (5o = corresponds to the 
constant /q. This leads to five (real) parameters in the 
IR. With the eight parameters mentioned for the UV 
discussed earlier, we now deduce that the domain wall 
solutions will be specified by a single parameter which 
we can take to be the wave-number k. 

We have constructed these solutions for the range 
0.253 < k < 0.75, which required utilising high precision 
numerics. As noted above we find that the T = limit of 
the black hole solutions approach these domain wall solu- 
tions. For example, in fig. 1 the blue line depicts the free 
energy of the domain wall solutions for various values of 
k, showing precise agreement with the T = limit of the 
black holes. Similarly, in fig. 2, for the domain wall with 
/c = fco ~ 0.256 we have shown the UV values of c„ and k 
with a blue dot and again we see precise agreement with 
the corresponding black hole solution. 



FINAL COMMENTS 

Let us summarise the main physical results. The d — 4 
CFT, dual to the AdS^ solution of our model (IT]), held at 
finite chemical potential undergoes a second order phase 
transition at a critical temperature Tc. The new phase 
is a helical superconducting phase that spontaneously 
breaks both the global abelian symmetry and the three- 
dimensional spatial Euclidean symmetry down to Bianchi 
VIIq symmetry. At Tc the spatial modulation is fixed 
by a wave number fc^ and as the temperature is cooled, 
the wave-number monotonically decreases. The T = 
ground state of the system maintains the helical order 
with non-vanishing wave number and has an emergent 
scaling symmetry in the far IR. These homogeneous, non- 
isotropic ground states at T = are holographically de- 
scribed by smooth domain wall solutions. It is natural to 
next use the AdS/CFT correspondence to analyse trans- 
port. Based on the rich optical properties of other helical 
orders (such as the chiral nematic phase of liquid crystals, 
recently discussed in |14j ) we expect interesting results. 

We expect analogous solutions for the helical supercon- 
ductors of the D = 5 models with SU{2) x t/(l) gauge 
fields studied in [9]. In particular we have checked that 
the (top-down) Romans' theory admits a Bianchi VIIq 
ground state solution with scaling exponent z w 3.98 
which we conjecture to be the IR limit of a T = do- 
main wall solution. More generally we now anticipate 
many other constructions of black holes in D = 5, 4 that 
are dual to other spatially modulated phases, both su- 
perconducting and otherwise. 
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